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We express the low-frequency shot noise in a disordered normal-metal - superconductor (NS) 
junction at finite (subgap) voltage in terms of the normal scattering amplitudes and the Andreev 
reflection amplitude. In the multichannel limit, the conductance exhibits resonances which are 
accompanied by an enhancement of the (differential) shot noise. In the study of multichannel single 
and double barrier junctions we discuss the noise properties of coherent transport at low versus high 
voltage with respect to the Andreev level spacing. 



Shot noise arises from the current fluctuations in trans- 
port as a consequence of the discrete nature of the car- 
riers, first predicted by Schottky for the vaccuum tubeEI. 
In a coherent conductor the quantum fluctuations, which 
follow from the probabilistic nature of the backscattering 
restricted by the Pauli principle, are reduced in compar- 
ison to the Schottky resultn. In a disordered normal- 
metal - superconductor (NS) junction, the shot noise is 
produced by the normal scattering-processes as well as 
the imperfect Andreev reflectionSij. The Andreev re- 
flection introduces fluctuations proportional to the dou- 
ble electron charge, which may be iteratively increased 
to give fluctuations of several charge quanta in biased 
SNS junctional In this paper we study the shot noise 
at finite voltage in a disordered NS junction as well as 
its structure due to the iterative scattering processes be- 
tween the disordered normal lead and the NS interface. 
We derive a general formula which expresses the differ- 
ential shot noise of the dirty NS junction in terms of the 
scattering matrix of the normal lead and the (scalar) am- 
plitude of Andreev reflection. In the models of a single 
barrier normal-metal - insulator - normal-metal - super- 
conductor (NINS) junction and a double barrier NINIS 
junction, we explain the existence of a non-trivial reso- 
nance structure. 

We consider a disordered NXS junction, shown in the 
inset of Fig. [|, with an arbitrary elastic scattering re- 
gion X in the normal lead, whose effect on the noise is 
to be determined. The low-frequency power spectrum of 
the current fluctuations is determined by the irreducible 
current-current correlator 



can be expressed through the positive energy eigenfunc- 
tions, 



P= / dte^((I(t)I(0))), 







(1) 



(((/ (t) I (0))) = ((/ (t) - (/)) (/ (0) - (J))) is the second 
cumulant in time) . The time-dependence of the current 
operator is defined by 

I (<) = exp [i (H - fiN) t] I exp [-i (H - fiN) t] . (2) 

In the mean field approximation for the Hamiltonian H 
we do not account for the fluctuations of the order pa- 
rameter in the S region. The effective Hamiltonian is 
diagonalized by a Bogoliubov transformationtl. After in- 
tegration over the cross-section, the net current operator 



I(t) 
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which may be evaluated in the normal lead. The oper- 
ators 7„j belonging to the wavefunctions {u m ,v m ) anni- 
hilate the scattering states of the disordered NS junction 
(note ud x v = ud x v — vd x u). We note that the symme- 
try of the BdG eigenfunctions (u n ,v n ), e„ with respect 
to sign reversal of energy to (— u* , u* ), — e n is the conse- 
quence of spin degeneracy in Nambu space, allowing the 
the negative energy states to be eliminated in favor of 
the positive energy states. 
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FIG. 1. Differential conductance (top) and shot noise (bot- 
tom) in a multichannel NINIS junction (see inset, symmetri- 
cal barriers of strength J dxV (x) — ZUvf, Ti = T2 « 0.05, 
interbarrier distance d = 20«f/A, 8 x 10 4 channels). The 
average shot noise per channel [Eq. (^)] exhibits maxima at 
the resonances of the conductance (solid lines). Note the en- 
hanced structure of the noise with respect to the conductance 
[Eq. (g)]. The dashed lines show the corresponding results 
for the normal state NININ junction. Inset: schematic NXS 
junction, e.g., X = INI. 

At voltages below the superconducting gap (eV < A), 
the quasi-particles are injected from the normal reservoir 
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only, and the wavefunctions depend on the N x N re- 
flection matrices r ee (e), r he (eX, r eh (e), and r hh (e) of the 
entire disordered NS junctional. The scattering states in- 
dexed bym = [|j,/i,6] consist of an incident e(lectron) or 
h(ole) like quasi-particle in channel /x with energy e su- 
perposed with the reflected electron and hole states. The 
occupation numbers are given by the Fermi-Dirac distri- 
bution f e h — f (e ± eV) (the voltage is measured with 
respect to the chemical potential of the superconductor) . 

The noise power (Q) is determined by the transitions 
induced by the current operator / from initial states 
\n) = \n n = l,n m = 0) to intermediate states |m) — 
|n„ = 0, n m = 1) and back, which differ only by their oc- 



cupation with respect to two single particle states with 
indices n and m. E.g., the transition between an in- 
cident electron n = [e, v, e] and an incident hole state 
to = [h, /i, e] is produced by the matrix element between 
the reflected electrons and holes of the two scattering 
states with respect to the current operator, (77i|/|?i) oc 

fe (1 - fh) (r eh r e e - Thh r he)iiv- Summed over the chan- 
nels, the transitions contribute to the fluctuations with 
the weight E^IH^)! 2 oc / e (l-A)Tr{(rt e r e/l - 
rl e r h h){r\ h r e e ~ r{ h r he )} = f e (1 - f h )Tr{rl e r he (l - 
r jj e r he)}- Following the quantum mechanical formalism 
outlined, rwe obtain the low-frequency limit of the power 
spectrumliSI valid for T <C eV < A, 
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[fe (1 - fh) + fh (1 - fe)] ^[rljhe (l - r{ e r he ) ] + [f e (1 - f e ) + f h (l - fh)] Tr [ (r{ e r he ) 



(4) 



The first term describes the transitions between states 
of opposite current sign, while the second represents the 
transitions between states of equal current sign. At zero 
temperature, the second term is suppressed by the PAi4i 
principle and the first term produces the shot noiscBO. 
The thermal fluctuations describing the Johnson- Nyquist 
noisetj are due to both terms. Eq. ( 41 is manifestly in- 
variant under sign reversal of voltageEl 

Interplay between normal scattering and Andreev re- 
flection. We express the electron-hole reflection matrix 
r/j e (e) in Eq. (^) in terms of the normal reflection- and 
transmission matrices ru and tij of the scattering (X) re- 
gion and the amplitude of Andreev reflection F at the XS 
interface, r he (e) = t* 12 (-e) [l - T 2 (e) r 22 (e) r* 22 (-e)] ~ 
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V ( £ ) = ( r L r he) 



r (e) t 2 \ (— e) (scattering matrices for X, see Fig. [l]: 
ru(e), t^ (e) for electrons, r*j (— e), t*j (— e) for holes, 

i,j = 1(2) left (right), r(e) = e M ^\ d{e) = 
arccos (e/A)). Inserting r^ e (e) in Eq. |] provides us with 
a general multichannel expression for the shot noise in a 
disordered NXS junction, explicit in the normal scatter- 
ing matrix and the Andreev reflection amplitude. 

Next we restrict ourselves to the situation of a junc- 
tion with uniform transverse cross-section, which permits 
the separation of channels and will be used inJhe single 
and double barrier systems below. We recalls that the 
differential conductances of the channels v — 1...N are 
given by the matrix product r he The in Eq. (^J), 



T v {e)T v {-e) 



T v (e) T v (-e) + R v (e) + R v (-e) - 2Re T (e) 2 r v (e) r* (-e) 



(5) 



r„ (e) is the reflection amplitude of channel v in the ma- 
trix T22 (e) describing the reflection from the right-hand 
side of the X region, R u (e) = \r v (e) | 2 , and T v (e) = 
1 — R v (e). The denominator of Eq. (g) contains the 
interference term describing the iterative scattering pro- 



cesses between the X region and the XS interface. At 
zero temperature, the shot noise (|j) exhibits the noise 

power P = - Jq'^' de ^„ S v (e) with the differential 
(low-frequency) shot noise of channel v, 



^S v (e) = 4 {r{ e r he [l - r{ e r he 
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The two energy dependencies of the reflection proba- 
bilities R v (e) and (the phases of) the reflection ampli- 
tudes r„ (e) translate into the voltage dependence of the 
shot noise. The energy dependence of the Andreev am- 
plitude can often be neglected. In the limit e — » 
(r — > — i) the phase dependencies drop-out of Eq. (||) and 
we recover the linear response result^ (h/2e 3 ) S v (0) = 
16T(0) 2 [l-T(0)]/[2-T(0)] 4 . 



NINS junction. In order to decribe the weak cou- 
pling of a normal lead to an NS proximity sandwich, 
we consider a normal-metal - insulator - normal-metal 
- superconducting (NINS) junction with a potential bar- 
rier of low transmission (T < 1) placed at a distance 
d away from the perfect NS interface. The NINS junc- 
tion serves as a model system for a tunneling experi- 
ment from a metallic tip to a thin film NS layer struc- 
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ture, which has permitted the observation of the Rowell- 
McMillan oscillationsEj. In this elementary model, the 
reflection amplitudes r v (e) = \fRe lVv ^ have a roughly 
constant modulus R and an energy dependent phase 
tp v (e) ~ 2 (k v + e/v v ) d + ipo, accumulated during the 
propagation between the potential barrier and the NS 
interface (k u , v v denote the Fermi wave number and ve- 
locity of channel v) . Consequently the shot noise (^|) for 
the channel v , 



2e 6 



4T 2 2i?[l - cosa„ (e)] 
(T 2 + 2i?[l-cosa„ (e)]) 2 



(7) 



depends only on the phase difference a v (e) — <p v (e) — 
(p v (— e) — 2$ (e) of electron and hole propagation through 
the X (=IN) region. The resonance structure of S u is 
intimately connected to the voltage (energy e = eV) de- 
pendence of the differential conductance (||) , 



2R\\-cosa v (e)]' 



(8) 



Figures ||a) an d @(b) show their generic dependence on 
a v . The minima of the denominator [cosa„ = 1 => 
(-ij.n = v v /2d(mr + arccos e„ )n /A)] correspond to the res- 
onances of the Andreev quasi-bound states. These con- 
ductance resonances are repelled from zero voltage, since 
the phase of the Andreev reflection d (e — > 0) = tt/2 has 
to be compensated by the phase difference of electron 
and hole propagation tp v (e) — (p v (— e). The shot noise 
(0) vanishes at these resonances (£„ {e u<n ) = 0). The 
noise maxima are found from cosa„ = (2R — T 2 ) /2R 
at energies doubly peaked close to the resonances. The 
peak separation Se ~ Tvp/d coincides with the width of 
the resonances. 

Interestingly, this non-trivial structure survives in the 
multichannel NINS junction, whose results are shown in 
Figs. ||(d) and|^(e). The stability of the resonance struc- 
ture is due to its pinning to the Fermi energy. In contrast, 
the resonance structure of a potential well in a normal 
conductor is washed out in the multichannel limit. As 
Figures ||(b,e) show, the narrow double peak structure in 
the noise is smeared out in the multichannel limit, and 
the noise S takes a maximum rather than a minimum at 
the resonance, retracing the shape of the conductance. 

At large voltages (eV ^> Vp/d), the shot noise ap- 
proaches a constant as a consequence of the dephasing 
between the channels. The channel average may be eval- 
uated by averaging over the phase a v in Eqs. (0) and (||) 
and we obtain the result (a„ (l/AT^ n — > l/2n da), 



S=±- 
2tt 



27T 2e 3 - 2e 2 

da S v (a) = —T, G = —T. (9) 



It follows that both the shot noise and the conductance 
approach the normal state values at voltages vp/d <C 
eV -C A. This limiting behavior demonstrates that the 



NIN junction is effectively decoupled from the NS inter- 
face in the large voltage limit and dominates both noise 
and conductance due to its low transparency. 
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FIG. 2. Zero temperature conductance (a,d), differen- 
tial shot noise (b,e), and finite voltage noise power to 
current ratio (c,f) for a NINS junction with a barrier of 
strength J dxV(x) = STlvf, average transmission T = 0.05, 
<lf = 500A, d = 20«f/A (v F /d < A < e F ). (a,b,c) fol- 
low from Eqs. (Q) and (|J) for one channel. (d,e,f) from Eqs. 
(^) and (^|) are averaged over 8 x 10 4 channels. Note the 
pronounced resonance structure in both conductance (d) and 
noise (e), which approaches the values for the correspond- 
ing NIN junction at large voltages (dashed lines). The noise 
to current ratio (c,f) decays to the classical Schottky value 
P/el = 1 above the Andreev minigap. 

Let us concentrate on the overall ratio of noise power 

to current P/el [P = f^ Vl deS{e), I = j* V deG(e)}. 
This ratio has already been used successfully in de- 
termining the unit of the charge carriers in the Frac- 
tional Quantum Hall Effectlij. While at small voltages 
(P/el = 2) the noise carries the signature of the Cooper 
pairs created by the Andreev reflection, at large voltages 
{vp/d <C eV -C A) it decays to the normal state value 
(P/el = 1). Interestingly, this decay is immediate at 
the first Andreev resonance, as is seen from Figs. 2(c) 
and (f). In fact, in a single channel, the energy inte- 
gration is equivalent to the integration over the phase 
a„, which at the midpoint between two resonances a v 
has covered the period 27r, thus producing P/el = 1 at 
a v = 7r, see Fig. 2(c). In the transport through a normal 
lead weakly coupled to a NS sandwich, we thus find two 
characteristic regimes: At voltages below eV < vp/d the 
one-particle excitations in the NS structure are effectively 
gapped and we may only couple to the superconducting 
condensate which introduces the charge 2e in the noise 
power. At voltages above the first Andreev resonance, 
we may couple directly to a one-particle density of states 
in the normal layer and we loose the signature of the 
superconductor at once. 
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Double barrier junction. In the study of normal versus 
superconducting junctions, it is of interest to investigate 
the impact the quantum coherence has on the macro- 
scopic limit of transport (limit of an infinite number of 
channels) and, in particular, to compare the behavior of 
a specific scattering region connected to cither a normal 
or a superconducting reservoir. We consider an insula- 
tor - normal-metal - insulator (X=INI) interferometer 
which we connect to normal leads in a NININ junction 
and to a superconducting lead in a NINIS junction. Be- 
sides the intrinsic interest in the double barrier system as 
a paradigm, now to be combined with a superconductor, 
it may also serve as a qualitative model for a disordered 
conductor, due to its strong similarity in the transmission 
distributionO. 

The resonance structure of the INI interferometer re- 
sults in a bimodal distribution of its (overall) transmis- 
sion probabilities T. For symmetric barriers it takes the 
form (infinite channel limit, T\ — T 2 <C 1)lj, 



P(T) = 



lTi 



1 



7T 2 y?T 3 (1 - T) ' 



re 



n 



(10) 



which has its analog in a disordered conductorlizl. We 
recall that, in a NININ junction, the rich structure 
of this distribution has no impact on the macroscopic 
transport properties. The (linear) conductance and shot 
noise may be evaluated from the bimodal distribution 
©, {h/2e 2 )G = JdTp(T)T = Ti/2, (h/2e 3 ) S = 
JdTp(T)T(l-T) and show no characteristics 

of coherent transport^. The conductance follows from 
the series resistance of the two barriers, and an incoher- 
ent model of a double barrier junction gives the same 
suppressed noise S/eG = 1/2 as a consequence of charge 
conservation^. When considering the INI interferometer 
in a NINIS junction instead, we find a non-trivial answer, 
both at small and large voltages. The linear response fol- 
lows from Eqs. (g), and (10), 
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27 
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~2e* 



2 G(0) 



S(0) 



dTp{T) 



dTp{T) 



2T 2 
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16T 2 (1 



2 V2' 
T) T x 



(2 - T) 



2 4V2" 
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The ratio of shot noise to conductance isfca 
S(0)/eG(0) = 3/4. A comparison of the shot noise 
to conductance ratios in a single barrier junction (X=I 
or X=IN) and a junction with disorder (X=D) is in- 
structive, see Table 1. In the single barrier junction 
the distribution of transmissions is peaked at T <C 1. 
The noise is due to the Schottky type fluctuations in 
nearly closed channels, and consequently the noise ratios 
S N /eG N = 1 in the NIN and S s /eG s = 2 in the NIS 
junction differ only by the charge quanta involved. In a 
junction with disorder, the noise-to-conductance ratio is 
also doubled for the NDS case, see Table 1. However, 
since in the bimodal distribution of the transmissions 



the noise is produced by the channels with intermediate 
transmission < T < 1 and the current by the open 
channels with T — > 1, the doubling is a non-generic and 
thus yet unexplained feature. This is demonstrated by 
the noise-to-conductance ratio of 3/4 in the NINIS as 
compared to 1/2 in the NININ junction. 

At finite voltage, the shot noise is described by Eq. 
(o) with the reflection amplitudes r v (±e) = r 2 + 
f |£ ie 2iM±e)d/ (l _ rir2€ 2iK(±z)d^ of the double barrier 

interferometer. We have evaluated this expression for 
a symmetric double barrier junction (Ti = T 2 = 0.05) 
and display the results in Fig. [[} At voltages of the or- 
der of the Andreev levels eV ~ Vir/d, we find a reso- 
nance structure independent of the number of channels. 
We observe again that the differential noise follows the 
resonance peaks of the conductance. At large voltages 
(eV ^> vp/d), the resonances disappear as the electrons 
and holes dephase, and we expect the approach to a 
regime where the conductance and shot noise become in- 
distinguishable from the incoherent addition of the NIN 
and NIS junctions. 

In conclusion, we have expressed the differential shot 
noise in a disordered NXS junction in Eq. (||) in terms 
of the normal and Andreev scattering amplitudes. We 
have described the resonance structure in the shot noise 
found at finite voltage in coherent transport. The ro- 
bustness of the resonance structure in the multichannel 
limit is owed to its pinning to the Fermi energy. We have 
pointed out the possibility of a non-trivial normal ver- 
sus superconducting noise ratio in the NINIS junction as 
as consequence of the bimodal distribution. Finally, we 
have found a rapid decay of double shot noise in a NINS 
junction above the Andreev minigap. 

We acknowledge fruitful discussions with D. Agter- 
berg, Ya. M. Blanter, M. Dodgson, and V. B. Geshken- 
bein. This work has been supported by the Swiss Na- 
tional Science Foundation. G. B. L. acknowledges partial 
support by RFFRN 96-02-19568. 

TABLE I. Shot noise to conductance ratio S/eG of a NXN 
compared to a NXS junction; single barrier X = I, double bar- 
rier X = INI, disorder X = D. The results are validpfotiStnaJl 
transmissions T <C 1, and many channels (N — > co)HoEj , Lj. 





single 


double barrier 


disorder 




barrier T < 1 


Ti = T 2 < 1 




S'jv/eGjv 


1 


1/2 


1/3 


Ss/eGs 


2 


3/4 


2/3 
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APPENDIX: 

We diagonalize the BCS-Hamiltonian for supercon- 
ductivity by a Bogoliubov transformation, considering a 
spatially inhomogeneous system. The transformation is 
carried out by determining the quasiparticle wavefunc- 
tions that fulfill the Bogoliubov-de-Gennes (BdG) equa- 
tions. We show that the symmetry in the solutions of 
the (BdG) equations with respect to reversing the ener- 
gies e a — ► — e a is a consequence of the spin degeneracy. 
Making use of the spin-reversal symmetry, we express 
the Hamiltonian, the density, and the current operators 
in terms of the quasiparticle operators 7 Q , and arrive at 
the current expression used in Eq. (^) . 

The mean- field, spin singlet BCS Hamiltonian for su- 
perconductivity can be expressed in the form (ho = 
(-iV + eA) 2 /2m) 

H = jd* X :& (x) ( ho ~ * _ h f + ^ ) * (x) : (Al) 
using the Nambu spin-up annihilation operator 



*(x) 



(A2) 



The pair potential is given by A(x) = X{^\ (x) (x)) 
(coupling constant A), the colon (:) denoting normal or- 
dering with respect to and The Hamiltonian ( Al) 
can be diagonalized by a basis transformation, 



*(*) =^*«(x) 7ct , 

a 

7a =y d 3 x*t ( x )#(x), 



(A3) 
(A4) 



with the eigenfunctions $ Q (x) of the BdG equations 
(which follow from the insertion of the Hamiltonian ( Al) 
into [H, 7«] 



h - n 
A* 



A 

-h% + n 



$ Q (x) = e Q * Q (x) , 



(A5) 



*a (x) 



Note that j a are annihilation operators of spin-up states 
and thus the excitations are all described in terms of spin- 
up quasiparticles. In order to preserve the (fermionic) 
commutation relations we need a complete orthonormal 
set of wavefunctions of the hermitian operator in (A5), 



£* a (x)* a (0 =15(x-x') 

a 

J d 3 x (x) ® a , (x) = 5 a , a , 



(A6) 



(A7) 



(involving n both positive and negative energy 
eigenstates)tl. 

Consider the spin-reversal transformation S, 



S^S 1 = 





-I 



(*.)■ 



(A8) 



S linear (Stp^S^ 1 —ipi, Si/j^S^ 1 = — V'f)- Noting that the 
order par ameter A(x) is invariant under the transforma- 
tion (A8), it is easily seen that the Hamiltonian (Al) is 
spin-reversal symmetric, [H, S] = (this symmetry ex- 
tends also to finite magnetic field, if the Zeeman splitting 
is neglected) . By means of spin- reversal we may attribute 
to each quasi-particle operator 7 a a linearly independent 
operator "f & through 



(A9) 



4 = s^s- 1 = J d\ *t (x) (o -i^ K (x) 

7a describes an excitation with opposite energy 
(according to [H,ja] — £c*7a)- From Eqs. (A9) and ( A4) 
we infer the effect of spin-reversal on the electron-hole 
wavefunction, 



-1 

1 



(A10) 



We arrive at a complete set of quasi-particle excitations 
in spin- up space, which are grouped into pairs j a , 7 a with 
energies ±e Q , as a direct consequence of spin degeneracy. 
Note that within the Nambu picture, all quasiparticles 
carry spin up [spin f e and spin f h (instead of spin J. e)], 
explaining the opposite energy of the related wavefunc- 
tion ( Alpp . The Hamiltonian (Af ) takes the form, 



n= J2 e - 



e a >0 



ilia 



d 3 x \v a \ : 



(All) 



The ground state is realized by filling all the negative 
energy quasi-particle states. 

The spin-reversal symmetry allows us to express all 
equations using only half of the eigenstates, i.e., one rep- 
resentative from each pair of states {$> a , 3?a}- In the 
following we choose the positive energy eigenstat es, ex - 
pressing the negative energy eigenstates through ( |A1C| ). 
We keep the positive energy states for the spin-up quasi- 
particles ("f n = 7n|), and reinterpret the related quasi- 
particle states of opposite (negative) energy as spin-down 
excitations (jn = 7^ ,). The Bogoliubov transformation 
then takes the well-known form 



7«T 



(A12) 



and the Hamilton operator is expressed by TL = 

T,<,, 1>0 e n(7„ T 7nT +7iu7nJ. ~ 2 / d 3 x \v n \ 2 ), displaying the 
spin degeneracy in the usual fashion. In the same way, 
we give the density and current operators in both re- 
spresentations, using all (indexed by a) or only the 
positive energy eigenstates (indexed by n), respectively 

(u V v = uS/v — (Vu)v), 
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p(x) =-e:*t( x) M ^ j*( x ):=- e ^«(xK,(x) 7 t 7Q , +u ;(xK,(x)7a'7i) (13) 
= -e ^2 [ w m( x ) M »( x ) ~ u m( x H( x )] (7™ T 7«T + 7^;7n|) ~ 2e l w "( x )| 2 

e m ,e„>0 e„>0 

-e ^ {[w m (x)w n (x) + u m (x)u„(x)] 7 m t7fii + h.c. } , (14) 

e m ,e„>0 



j( x ) 



T : *t (x) V* (x) 



e 
2mi 



J2 u * a ( x N u a'i x hha' -<(x)V«a'( X )7a'7l (15) 
a, a' 

1 H ( M m( x )Vl»W + ?4(x)V«n(x)J (7„ T 7nT + 7lx7«i) ~ 2 ^ <( X )V^n( x ) \ 
-. ^ |(^ M m(x)V!ln(x) + ll m (x)\7« n (x)J 7 m |7s 1 - Il.C. j. (16) 



e 
2mi 



e m ,e„>0 



The current operator is easily generalized to finite mag- 
netic field. We note that current and density operators 
obey the continuity equation 



0p(x) 
dt 



V • j(x) = 2ie* f (x) 







A(x) 
-A*(x) 



*(x), 
(17) 

where the right hand side vanishes when taking the (self- 
consistent) expectation value. The effective current ex- 
pression used as a starting point in Eq. (||) is obtained 
from (|l6|), accounting for the spin degeneracy by a fac- 
tor 2 and dropping the particle non-conserving last term, 
which is eliminated by the frequency integration in ([!]). 
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